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Abstract 

We present a unifying treatment of dark energy and modified gravity that allows 
distinct conformal-disformal couplings of matter species to the gravitational sector. 
In this very general approach, we derive the conditions to avoid ghost and gradient 
instabilities. We compute the equations of motion for background quantities and linear 
perturbations. We illustrate our formalism with two simple scenarios, where either 
cold dark matter or a relativistic fluid is nonminimally coupled. This extends previous 
studies of coupled dark energy to a much broader spectrum of gravitational theories. 
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1 Introduction 


The nature of dark energy, responsible for the present cosmological acceleration, is a 
central topic in theoretical and observational cosmology. One of the main goals of current 
and future cosmic surveys is to constrain or possibly detect deviations from the standard 
ACDM scenario, induced by the presence of dark energy or modifications of General 
Relativity (GR) (see e.g. [1]). This is particularly relevant on scales above ~ 10 Mpc, 
where deviations from GR are not yet well tested. Fortunately, on these scales cosmological 
perturbations are still in the linear regime today and linear perturbation theory around a 
FLRW background is thus a valid description. 

Given the plethora of existing dark energy and modified gravity models (see for in¬ 
stance [2,3]), it is worth resorting to an effective approach that tries to describe all possible 
deviations from ACDM in a simple and systematic way, relying on a minimal number of 
parameters. In the linear regime for perturbations, this task has been sucessfully under¬ 
taken for single scalar field models in [4,5]. Initially inspired by the so-called Effective 
Field Theory of inflation [6,7] and minimally coupled dark energy [8], this approach relies 
on the construction of an effective action for linear perturbations. In order to do so, we 
start from a generic Lagrangian written in terms of Arnowitt-Deser-Misner (ADM) [9] 
quantities defined with respect to the uniform scalar field hypersurfaces (see also [10,11] 
for an analogous approach, [12-14] for recent reviews and e.g. [15-17] for applications). 
After having been implemented in a public numerical code named EFTCAMB [18], most 
recently, it has been applied to constrain deviations from the standard cosmological con¬ 
stant scenario by the Planck collaboration [19]. 

The action developed in [5] contains five free functions of time that parametrize any 
deviation from ACDM. Four of these functions describe cosmological perturbations in 
Horndeski theories [20-22], The fifth parameter describes deviations from GR encompass¬ 
ing Horndeski theories. Indeed, the same formalism was also instrumental to uncover the 
theories beyond Horndeski of [23,24], which lead to equations of motion higher than sec¬ 
ond order but are free from Ostrogradski instabilities (see e.g. [25] for an earlier example 
of theories beyond Horndeski). 

The developments described above assume that matter is minimally coupled to a unique 
metric, which will be called Jordan frame metric for convenience. However, although the 
universality of couplings is very well tested on Solar System scales [26], on cosmological 
scales constraints are much weaker. In particular, the scalar field responsible for the 
current accelerated expansion is known to mediate a fifth force [27], which may lead to 
violations of the equivalence principle (EP) on large scales [28] (see also [29] for a test of 
the EP on large scales). Moreover, while fifth force effects on standard matter such as 
baryons and photons are severely restricted, those on cold dark matter (CDM) or neutrinos 
could be much larger. This leaves the freedom to consider the case where different matter 
species 1 couple differently to the scalar field [30,31]. 

The goal of the present work is to extend the approach developed in Refs. [4, 5] by 
relaxing the assumption that all matter species are minimally coupled to the same metric. 2 

1 By matter species we intend the different components in the Universe (baryons, photons, CDM and 
neutrinos) but the results derived here could be straightforwardly extended to different types of objects, 
such as e.g. galaxies of different sizes, behaving differently under the effect of the fifth force. 

2 Another general parametrisation of theories of single-field dark energy that is explicitly coupled to dark 
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For simplicity, in the following we restrict our study to effective theories of dark energy 
or modified gravity that remain within the Horndeski class. This means that we assume 
that the function oh introduced in [24] to describe theories beyond Horndeski at the level 
of linear perturbations vanishes here, leaving only four out of the five free independent 
functions of [5]. We reserve a treatment of theories beyond Horndeski for future work. As 
shown in [35], the structure of the Horndeski Lagrangians is preserved under a disformal 
transformation [36] of the metric with coefficients that depend only on the scalar field (not 
on its gradient), i.e. of the form 

Qtiv = C{4>)g^v + D(<j>)d il 4>d v ^ . ( 1 . 1 ) 

Thus, in the following we assume that each matter species is minimally coupled to a 
distinct Jordan metric of this form. * 3 While conformal couplings (i.e. with D = 0) have been 
extensively studied in the literature (see e.g. [1] and references therein), disformal couplings 
have been investigated only recently (see also e.g. [39-48]). Moreover, the dynamics of the 
gravitational metric g^ u is usually assumed to be governed by the standard Einstein-Hilbert 
action. Here, we allow a much more general gravitational sector, based on the effective 
description given in [5]. In Sec. 2 we review our formalism within the ADM effective 
approach and the gravitational action in the uniform scalar-field gauge. Apart from the 
four time-dependent parameters mentioned above, we introduce two extra functions of 
time for each species, describing the nonminimal coupling to dark energy via an effective 
metric of the form (1.1). The structure of this action is preserved under transformations 
of the reference metric of the form (1.1) and the stability conditions for the matter and the 
gravitational sector are shown to be invariant under these transformations. More details 
on the frame dependence and on the derivation of the stability conditions of gravitational 
and matter quantities are respectively given in Appendix A and Appendix B. 

In Sec. 3 we derive the evolution equations describing the matter sector, which now in¬ 
clude the effect of the nonminimal couplings, and in Appendix E we provide the definitions 
of several parameters introduced in this section. These equations must be supplemented 
with the Einstein equations describing the gravitational sector, reported in Appendix C. 
We provide and discuss the perturbation equations using Newtonian gauge but these are 
also given in synchronous gauge in Appendix D. 

The parameters of our effective description can be constrained by observations. As a 
direct application of our approach, in Sec. 4 we consider the cosmological consequences, for 
the background evolution and for linear perturbations, of a Universe where the coupling 
of CDM differs from that of the other species (see e.g. [49-59]). Our analysis extends 
previous results as we allow gravity itself to be modified, not only the couplings to matter. 
In Sec. 5 we consider the case where the coupled species is a relativistic fluid. This will 
allow us to highlight the dependence of conformal and disformal couplings on the equation 
of state. Finally, we conclude in Sec. 6. 

matter has been given in [32], in the framework of the Parameterized Post-Friedmann approach [33,34]. 

3 Other types of couplings can be found in the literature. For instance, Ref. [37] considers a CDM 
action that depends on the contraction of the CDM 4-velocity with the normalized space-time gradient 
of the scalar field, in the context of Lorentz-violating theories. Ref. [38] directly modifies the action for a 
general perfect fluid. 
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2 Unifying description of dark energy with non universal 
couplings 

In this section we introduce the gravitational and matter actions within the ADM frame¬ 
work introduced in [5] and recently summarized in [14]. After giving the background 
equations of motion, we study linear fluctuations and derive the conditions for the linear 
theory to be healthy, i.e. ghost-free and without gradient instabilities. 

2.1 Gravitational and matter actions 

In the present work, we assume that the gravitational sector is described by a four¬ 
dimensional metric g^ v and a scalar field (j>. Let us start by choosing a coordinate system 
such that the constant time hypersurfaces coincide with the uniform scalar field hypersur¬ 
faces. In this gauge, referred to as unitary gauge, the metric can be written in the ADM 
form, which reads 


ds 2 = —N 2 dt 2 + hij (dx l + N l dt) (dx- 7 + N^dt) , (2.1) 

where N is the lapse and N 1 the shift. In the following, a dot will stand for a time 
derivative with respect to t, and Di will denote the covariant derivative associated with 
the three-dimensional spatial metric hij. Spatial indices will be lowered and raised with 
the spatial metric hij or its inverse h lJ , respectively. 

In the unitary gauge, a generic gravitational action can be written in terms of geometric 
quantities that are invariant under spatial diffeomorphisms [6, 7]. Expressed in ADM 
coordinates introduced above, these geometric quantities are the lapse N, the extrinsic 
curvature K t j of the constant time hyper surfaces, whose components are given by 

Uj = 2jv(^b' — DiNj — DjNi ) , (2-2) 

as well as the 3d Ricci tensor Rij of the constant time hypersurfaces and, possibly, spatial 
derivatives of all these quantities. Thus, the action is generically of the form 

Sg = J d^xyf^gL(N,Kij,Rij,hij,Dijt) . (2.3) 

The gravitational action must be supplemented by a matter action S m . In order to 
describe dark energy and modified gravity scenarios with EP violations, we assume that 
beside the gravitational sector introduced above, the Universe is filled by Ns matter species 
labelled by an index I, with 1 = 1 ,..., Ns, each minimally coupled to a different metric. 
For each species I, we denote the corresponding metric by g^J and we call this the Jordan 
frame metric associated with this species. The total matter action is thus given by 

N S . _ 

S ra = ^S I , Sr = I d A x\J-gi 1 ) ^i) , (2.4) 

i J 

with 

9$ = (4>)9hu + D^ ( 0 )< 9 M 0 < 9 „<£ • ( 2 . 5 ) 
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In order to preserve the Lorentzian signature of the Jordan-frame metric of the species I, 
it is necessary to have > 0. 

There is some arbitrariness in the choice of the gravitational metric g since we work 
in the context of modified gravity, where the gravitational dynamics cannot be expressed 
in terms of a standard Einstein-Hilbert term, in general. It is often convenient to choose 
one particular matter species, say /*, and define its Jordan metric as the gravitational 
metric, in which case we have = 1 and = 0. 

2.2 Homogeneous equations 

Let us discuss briefly the evolution of the background metric described by a FLRW metric 
assumed to be spatially flat. In this case the lapse is a function of time only, which we 
denote N(t), the shift vanishes, N l = 0, and the spatial metric reads gtj = hij = a 2 (t)8ij 
where a represents the scale factor. Thus, the metric reads 

ds 2 = — N 2 (t)dt 2 + a 2 (t)dx. 2 . (2.6) 

The homogeneous dynamics depends on the gravitational Lagrangian L in eq. (2.3), 
which can be seen as a function L(N,a,a ) when the arguments are restricted to their 
background values, i.e. N = N, hij = a 2 (t)8ij, Rij = 0, and 

Kij = Kij = = Hhij , (2.7) 

where H = a/{aN) denotes the Hubble rate. Here and in the following, barred quantities 
are evaluated on the background. 

The variation of the matter action S m with respect to the metric g^ w defines the 
energy-momentum tensor, according to the standard expression 

V llJ = . ( 2 . 8 ) 

V~9 8g^v 

This definition applies even if the matter is minimally coupled with respect to a metric g^ v 
that differs from g jlv , as discussed in Appendix A. In the homogenous case, the energy- 
momentum tensor depends only on the energy density p m = —T ° 0 and the pressure p m = 
If there are several matter components, the previous quantities simply correspond, 
respectively, to the sums of the energy densities and pressures associated to each individual 
species, i.e. p m = Y, I Pi and p m = J2iPi- 

The background evolution equations are then obtained by taking the variation of the 
total homogeneous action 5 g + 5 m with the respect to N and a. As shown in [5], this leads 
to the equations 

L + NLjy — 3 HT = p m (2.9) 

and 

L-3HT-^ = -p m , (2.10) 

where the coefficient T is defined by 

<2 - n) 
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Equations (2.9) and (2.10) generalize the usual Friedmann equations. For GR, where the 
Lagrangian is given by L = Mp(K, L jK — K 2 + R)/ 2, one can check that the standard 
equations are recovered, since L = —3 MpH 2 , Ljv = 0 and F = —2 MpH. 

The Friedmann equations eqs. (2.9)-(2.10) can always be written as 

H 2 = (Pr n + PEE ) , (2.12) 

H + t^H 2 = (Pm +Pde) , (2.13) 

where M denotes the effective Planck mass, which can be in general time-dependent (it will 
be defined below from the second derivative of L with respect to the intrinsic curvature). 
The above equations can be interpreted as definitions of the homogeneous energy density 
and pressure of dark energy, respectively given by 

Pde = 3 M 2 H 2 — p m , pde = —M 2 (2H + 3 H 2 ) — p m . (2-14) 

These equations can also be shown to be equivalent to the Friedmann equations derived 
from the Lagrangian [4,5] 

L = — (4) R+ —-A, (2.15) 

where is the 4d Ricci scalar and c = c(i) and A = A (t) are time-dependent functions, 
respectively given by 

2c = /jde + Pde + H(M 2 ) — (M 2 ) , (2-16) 

2A = /jde — Pde + 5iL(M 2 ) + (M 2 ) . (2-17) 


2.3 Linear perturbations 

We now expand the gravitational action up to second order in perturbations, in terms of 
the perturbative quantities 

5N = N — N(t ), SKij = Kij - Hhij , (2.18) 


as well as Rij, which is already a perturbation since its background value vanishes. 

The second-order expansion of the gravitational Lagrangian involves first and second 
derivatives of L with respect to its arguments K^, Rij and N. It is convenient to introduce 
the time-dependent coefficients Q, Br, B , Ak, Ak, C, C, Ar and Ar respectively as 


dL 

dR^ 


G6] 


d 2 L 

dK 2 8K[ 
d 2 L 

dRj dR{. 
d 2 L 

dK{ dR’ k 


d 2 L „ d d 2 L 

dNdR) R i ’ dNdK) i ’ 

(2.19) 

■■Ak^S^ + Ak (Si 5$+ 5*6^ , 

(2.20) 

: Ar Sj Sf + Ar + , 

(2.21) 

-.C 8 ) 8 ^+cU\ 8 k j+ 6 ik 8 ^ , 

(2.22) 
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where all partial derivatives on the left hand sides are evaluated on the background. 
The form of the right hand side of these expressions is merely determined by the FLRW 
symmetries. The first and second derivatives of L with respect to the scalar N are simply 
denoted as Ljy and Tat at, respectively. 

In the following, for simplicity, we restrict our considerations to Lagrangians that lead 
to dynamical equations with at most two space derivatives. This is automatically ensured 
if we impose the conditions [5,14] * * * 4 

Ak + 2 Ak = 0 , C + C = 0 , 4 Ar + 3Ar = 0 . (2.23) 


We also impose the further condition 

Br = j*(A K ~ G ~ HC) , (2.24) 

which is equivalent to restricting the range of application of the expanded action to Horn- 
deski theories [5]. 5 

The second-order gravitational action can then be explicitly written in terms of all the 
coefficients introduced above. In fact, the quadratic action involves only a few combina¬ 
tions of these coefficients, which are represented by the following dimensionless parame¬ 
ters [14,17] 


_ 2NL n + N 2 L nn _ BN 

“ K =-2 WT k -’ “ B = Wm 


CTt = 


g + C/(2N) + HC 
Ak 


- 1 . 


(2.25) 


The effective Planck mass squared is defined by M 2 = 2 Ak- With this definition, M 
coincides with the time-dependent Planck mass introduced in eqs. (2.9) and (2.10) and in 
the action (2.15). Its possible time variation is characterized by 


1 d In M 2 
" M = NH dt 


(2.26) 


In terms of these parameters, one finds that the second-order gravitational action is given 

by 6 


= J d 3 xdt a 3 N^~ 


5K)5K{ - 5K 2 + R^r + (1 + a T ) 6 2 (VhR/c 


+ «K H 2 


N 

5 N \ 2 < ttctJ N 

) +4 a^HdR. 

N N 


(2.27) 


where d 2 denotes taking the expansion at second order in the perturbations. Moreover, we 
have omitted irrelevant terms that vanish when adding the matter action and imposing 
the background equations of motion. 

J Here we have corrected a typo in Ref. [14]. The coefficient in front of SK5R inside the bracket in 

eq. (55) (see v2 of the arXiv version) should be C/2, so that the condition in the second line of eq. (60) 

should read C* = C + C. With this correction, eq. (76) of Ref. [14] is equivalent to eq. (2.23) in this article. 

s To parametrize deviations from Horndeski theories at the linear level, Ref. [24] introduced the param¬ 

eter a H = (G + HC + NBr)/Ak — 1- Here we will assume qh = 0. 

6 To write this action, we have not assumed N = 1 as done in previous references [4,5,14]. In such a 
way the action remains explicitly invariant under a time reparameterization t —> t(t), which is convenient 
when changing frame. 









To verify that M plays the role of the Planck mass which canonically normalizes the 
graviton, let us write this action in terms of the tensor fluctuations, defined as the traceless 
and divergence-free fluctuations of the spatial metric, i.e. 


hij — o, (t) ( Sij + '7 ij) ; 'Yu — 0 — dm, ■ 
The above action then yields 



(2.28) 


(2.29) 


where the tensor sound speed squared is given by = 1 + «t- Absence of ghosts and 
gradient instabilities respectively require that the kinetic and spatial gradient terms are 
positive, i.e. that 

M 2 > 0 , a T > —1 , (2.30) 

which will be assumed in the following. 


2.4 Matter couplings and stability conditions 

To discuss the stability and determine the propagation speed of dark energy perturbations, 
one must also include quadratic terms that come from the matter action, because the latter 
depends on the gravitational degrees of freedom. In order to do so, we need to to take into 
account the fact that each matter species / is minimally coupled to a metric g^J defined 
in eq. (2.5). For later convenience, we define, for each matter species, the time-dependent 
quantity 

4> d\nCf ] 
acj = 2 HN d4> ’ 

which parameterizes how the conformal coupling affects physical observables; the impact 
of the disformal coupling is parameterized by the quantity 7 


oid,i 


(jt/NfDf* 

C^ ] - (0/ N) 2 Df ] ’ 


(2.32) 


and the right-hand side of these equations are to be evaluated on the background. Re¬ 
quiring that the Jordan frame metric is Lorentzian implies ao,/ > — 1 [35]. 

In unitary gauge, eq. (2.5) reads 


9 $ = Ci(t)g lx „ + D I (t)8°X, 


[I V 


with 


Ci{t) = > Dj(t) = <fi 2 (t)D ( f\<j ) (t)) . 

Then the parameters ac,i and ao,/ introduced above take the form 


«C,/ = 


1 d In Ci 

2 HN dt 


, «D,/ = 


Dr 


N 2 Cr - Dr ' 


(2.33) 

(2.34) 

(2.35) 


7 This parameter coincides with I/ 7 2 , where 7 is the so-called disformal scalar in the notation of 
Ref. [47]. 
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Combining the quadratic action for matter with eq. (2.27), one can extract a quadratic 
action that governs the dynamics of the gravitational scalar degree of freedom and the 
matter ones. The explicit calculation in the case of perfect fluids is presented in Ap¬ 
pendix B. The absence of ghosts is guaranteed by the positivity of the matrix in front of 
the kinetic terms. For the gravitational scalar degree of freedom, this condition is given 
by 


a = a K + 6a| + 3 ^ «d,/ flj > 0 , 

/ 

where we have introduced the (time-dependent) dimensionless density parameter 


(2.36) 


Ql ~ 3M 2 H 2 ’ 


(2.37) 


where we recall that M 2 is in general time dependent. As pointed out already in [43,60], 
the presence of a disformal coupling affects the ghost-free condition. 

For the matter sector, the analogous condition usually corresponds to the Null Energy 
Condition [61]. In the Jordan frame of each species /, this can be expressed in terms of the 
energy density and pressure by pi + pj > 0 (we use the symbol" to denote Jordan-frame 
quantities). In the frame of g^, this inequality becomes 


Pi + (1 + «D ,i)pi > 0, 


(2.38) 


where we have used that wj = (1 + a-£,j)wi (see Appendix A for various relations between 
quantities defined in distinct frames). 

The speed of sound for scalar perturbations can be read off from the quadratic action 
derived in Appendix B. One finds 


c s =- \ (l+a B ) 

a 


■^2 -aM+OT+aB(l+aT) 


+ 


[ 1 + ( 1 +«D j)wi 


0/ } , (2.39) 


where matter appears in the last term in the bracket, proportional to ^ ~2i(pi+Pi )• Absence 
of gradient instabilities is guaranteed provided that 

c 2 > 0 . (2.40) 


We also require that the propagation speed for each matter species, in its Jordan frame, 
is positive, c 2 1 > 0. 


2.5 Disformal transformations 

As mentioned earlier, there is some arbitrariness in the choice of the metric that 
describes the gravitational sector. Let us thus see how the description is modified when 
the reference metric undergoes a disformal transformation, of the form 

9ij,u 9fiu = C^Xfyg^ + . (2.41) 

In unitary gauge, this corresponds to the transformation 

9^u 9^v = Citygfiv + , (2.42) 
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with C(t) = CW> {(j){t)) and D(t ) = (0(f)) 0 2 (t). The effect of this transformation 

on the ADM quantities, on the background quantities and on the linear perturbations is 
described in detail in Appendix A. Here, we just present the main consequences on the 
parametrization of the couplings and of the linear perturbations. 

In analogy with (2.35), it is convenient to introduce the dimensionless time-dependent 
parameters 

aC = 2 HNC ’ “ D = N 2 C — D ’ ( 2 ‘ 43 ^ 


which characterize, respectively, the conformal and disformal parts of the above metric 
transformation. 8 

Let us first see how the gravitational action (2.27) changes under the transformation 
(2.42). As shown in Ref. [35], the structure of Horndeski Lagrangians is preserved un¬ 
der a disformal transformation. Indeed, using eqs. (A.2) and (A.3), one can check that 
(2.27) maintains the same structure with the time-dependent coefficients in the action 
transforming as 


M 2 


M 2 

C \/l + CUD 


(2.44) 


and 


a k 


a B 


a M 

At 


ax + 12as[ac + (1 + ao)^] — 6[ac + (1 + od)^] 2 + 30 m aD 
(1 + ac) 2 (l + «d ) 2 

1 + CUB ^ 

(1 + ac)(l + a©) 

«m — 2ac _ Ad 

1 + ac 2HN(1 + «d)(1 + «c) 

(1 + ot)(1 + «d) — 1 • 


(2.45) 


We can use these transformations, which depend on the two arbitrary functions ac and 
od, to set to zero any two of the parameters a a above. 

Finally, the conformal and disformal coefficients associated with the respective matter 
Jordan frame metrics are modified according to 


«D,7 

Ac,/ 


«D,/ — «D 
1 + «D 
«C,/ ~ «C 

1 + ac 


(2.46) 


Alternatively to setting two A a to zero, it is always possible to choose as the new reference 
metric g jiv one of the matter Jordan metrics, say g^c\ which then implies Ac,/* = Ad,/* = 

0. 

One can verify that all the stability conditions are frame independent. In particular, 
the quantities that appear in the no-ghost conditions, eqs. (2.36) and (2.38), transform as 


a 

(1 + ac) 2 (l + ao) 2 ’ 


Pi + ( 1 + A D ,/)p/ 


PI + (1 + «D ,l)P! 
C 2 (l + a D )V2 ’ 


8 We require C > 0 and a® > —1, see discussions respectively in Secs. 2.1 and 2.4. 


(2.47) 
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and since 1 + «d > 0 (see discussion in Sec. 2.4), their sign is indeed frame independent. 
It is also straightforward to check that all the propagation speeds, i.e. of tensor, scalar and 
matter fluctuations, transform in the same way and that their signs remain unchanged, 

= (1 + od)cji , (? s = (1 + OiT))Cg , (? s j = (1 + q;d)c^/ . (2.48) 

In summary, at the level of linear perturbations our gravitational sector is characterized 
by four time-dependent parameters «k 5 ctM and cep- Each species is characterized by 
two time-dependent parameters, associated with their conformal and disformal couplings 
respectively. A priori, for a system of Ns species coupled to different metrics, this gives 
a total of 2 Ns + 4 parameters. However, the general invariance of the system under 
an arbitrary change of frame, characterized by two parameters, reduces the number of 
independent parameters to 2{Ns + 1). 

In particular, action (2.27) can also be used to describe inflationary perturbations. 
In this case, matter can be ignored, i.e. N$ = 0, and one can always use eq. (2.44) to 
find a frame where the Planck mass is time-independent and cp = 1, without loss of 
generality [62]. Thus, inflationary fluctuations can be generically described in the frame 
where «m = 0 = cep by only two operators, those proportional to «k and ob, as in 
Refs. [6,7]. 9 


3 Matter equations of motion 


In this section, we leave the unitary gauge description introduced in the previous sec¬ 
tion, by “covariantizing” the action. This can be done explicitly by performing a time 
reparametrization of the form 

t <j> = t + ir(t, x) , (3-1) 


where the unitary time t becomes a four-dimensional scalar field </>. For convenience, we 
denote by n the fluctuation of </>. 

By substituting the above transformation into the total action S = S g + S m , we then 
obtain an action that depends on the scalar held (j) and an arbitrary metric g We 
will use this more general form for the action to derive the evolution equations for the 
gravitational and matter sectors. 

The equations of motion for the metric are obtained by varying the total action with 
respect to g^ u , 


5S 

$9^ v 


= 0 . 


(3.2) 


which provides the generalized Einstein equations. At linear order, they are explicitly 
given in Appendix C. 

To write the equations of motion for matter, we use the invariance of the matter action 
Sj under arbitrary diffeomorphisms, This implies 


v /t r (/) ^ + QjdA = o, (3.3) 

9 The sound speed of fluctuations in this case is (? s = — (2/a) [(1 + ob )(H/H 2 + cm) + oib/H] ; see 
eq. (2.39). Thus, for a constant a®, the usual gradient instability associated with the violation of the Null 
Energy Condition for H > 0 can be cured by requiring —1 < cm < —H/H 2 [6]. 
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where the function Qj, which characterizes the coupling between the matter species I and 
the scalar field, is defined by 

= ~ IN'W'w +v " • (3 ' 4) 

where a prime denotes a derivative with respect to </>. The expression on the right hand 
side is obtained by using the property that the matter action Si depends on the scalar 
field only through the Jordan metric eq. (2.5). 

Finally, the evolution equation for (j> can be obtained by variation of the total action 
with respect to 0, 5S/8cj) = 0. Thus, from eq. (3.4) we obtain 


1 8S g 




(3.5) 


In the following, we will study the above equations, first in the homogeneous limit and 
then restricting ourselves to their linearized version. 


3.1 Homogenous equations 


Let us first consider the homogeneous case, with the flat FLRW metric (2.6) where we set 
N = 1. The associated Friedmann equations are given in eqs. (2.9) and (2.10), or (2.12) 
and (2.13), with p m = Y,i Pi and Pm = J2iPi- 

For a FLRW background, the definition of Qj, eq. (3.4), reduces to 


Qi 


H PI 
1 + a d,j 


|«c,/ [1 — 3to/(l + an,/)] + «d,/ 



+ 277(1+ a D ,/) j ’ 


(3.6) 


where we recall that the conformal and disformal parameters ac,i and «d,/ are respectively 
defined in eq. (2.35). Substituting the above expression into eq. (3.3), one finds that the 
homogeneous matter evolution equation can be written in the form 


pi + 377(1 + wi - 7 /)/>/ = 0 , 


(3.7) 


where the dimensionless parameter 7 / is given by 

77 = x«C,/ [1 - 3ui/(l + c®,/)] - wianj + , D,/ -c . (3.8) 

3 QH{l + ar,,i) 

Taking into account (3.7), one can also check that Qi = 377 / 5 / 7 /. Note that the equation 
of state in the Jordan frame of the fluid I corresponds to wi = u>/(l + aD,/) (see Appendix 
A.l). Using this relation, one can check that for a relativistic fluid, i.e. wi = 1/3, the 
conformal term in (3.8) disappears, as expected from the tracelessness of its stress energy 
tensor. 

Given the Friedmann equations (2.12) and (2.13) as well as the continuity equation for 
matter, eq. (3.7), the homogenous energy density of dark energy satisfies 

Pde + [3(1 + wde) — «m] HpvE = H ^(om - 3y i)pi , (3.9) 

1 

where we have introduced the equation of state parameter for the dark energy component 
^de = Pde/pde- 
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3.2 Perturbation equations in Newtonian gauge 

We now consider a linearly perturbed FLRW metric in Newtonian gauge with only scalar 
perturbations, i.e., 


ds 2 = —(1 + 2 <&)dt 2 + a 2 (t)( 1 — 2'&)5ijdx l dx - 7 . (3.10) 

In this gauge, we decompose the scalar part of the stress-energy tensor for each species, 
at linear order, as 


T(l )°0 = ~(Pl + dpi) , 

(3.11) 

T (i)°i = Pd 1 + wi)diVi = -a 2 T (7) * 0 , 

(3.12) 

T (i) l j = ( Piwi + Spi)d l j + (d l dj - <ti , 

(3.13) 


where dpi and dpi are the energy density and pressure perturbations, vi is the 3-velocity 
potential and 07 is the anisotropic stress potential for the species /. As usual, we define the 
total matter quantities as dp m = J2i$Ph = J2i$Ph v m = Y)i(pi + Pi)vi/(p m + Pm) 
and cr m = J2i (J i- 

The continuity equation, for each species, can be derived from the time component of 
eq. (3.3). In Fourier space, at linear order, this reads 

k 2 

dpi + 3 H(dpi + dpi) - 3pi(l + wi)^ - pi( 1 + wi)—vi = Qitt + SQi , (3.14) 

a z 

where Qi and dQi are given respectively by eqs. (3.6) and (C.6). The space components 
of eq. (3.3) gives the Euler equation, which at linear order reads 

2 k 2 

pi{l+wi)vi+pi [wi - 3Hwi(l + wi)] vi+dpi+pi(l+wi)<&--—(Ti = —Qi [n+vi(l+wi)] . 

O CL 

(3.15) 

We can rewrite the equations above in terms of the density contrast <5/ = dpi/pi and 
using the explicit expression for dQi given in eq. (C.6). This yields 


fdv T \ k 2 

di + 3H(1 + a c ,i)(l + a D)I ) - widi - (1 + wi)-^vi = 3 [1 + (1 + a D .i)wi] T 

V Pi J a 

+ 2(1 + odj) [ac,/(l - 3'wj) - 3yj] - a D ,J ^ - vf + w i^ 2 7T 

- [2(1 + av,i)ac,iQ ~ 3wi) — 3wian,i — 37/(3 + 2aD,/)] Hi r 


+ 3 


ipiiH) + wianpH + ( ac,i + oid,/(1 + «c,/)) wiH 


7 T 


(3.16) 


and 


vi — 3H 


wi - 7 / 


wi 


vi T 


dpi 




k 5 


3H(l + wi)\ (1 + wi)pi 3(1 + wi)pi a 2 


oi = —3 H 


7/ 


-7T . 


1 + Wj 

(3.17) 
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As mentioned before, the equation of state parameter in the matter Jordan frame, wi, 
is different from the one in a generic frame, wj. This means that the relation between 
pressure and energy density perturbations depends on the frame. Indeed, because of 
the coupling to the scalar field, there is a non-adiabatic pressure perturbation [63] which 
appears in frames that are disformally distinct from the Jordan one (see also [64] for a 
similar remark). For an isentropic perfect fluid with c 2 „ j = wj = constant, this reads (see 
Appendix A.2) 


dpna,d,I = dpi 


Pig 

—dpi = PI 

PI 


2«d( < I > — 7r) + 


«D 

1 + a D 



(3.18) 


Let us comment on the initial conditions of the above equations. In the simplest 
case, one can assume that perturbations start in the adiabatic growing solution, which 
is justified if they have originated from single-field inflation (see e.g. [65]). In this case, 
their amplitude can be given in terms of the time-independent quantity 7£j n , defined as 
the long-wavelength limit (k -C aH ) of the total comoving curvature perturbation [66] 


TZ=-^ + H 


4> + H<5> 
H 


(3.19) 


In Ref. [14] it was shown that, in the absence of nonminimal couplings, the generalized 
Einstein equations and the evolution equations for the matter and field fluctuations admit 
the adiabatic solution 

= —(1 + (y.-r)TZ- m + (1 + «m )He —, fk = — TZi n + He , 

M 2 (3.20) 

dpi = -pie , dpi = -pie , vj = e , tt = -e , 

where 

e = —j a [M 2 ( 1 + ^ + <7 m ] dt . (3.21) 

One can check that these expressions are frame invariant and remain a solution even 
in the presence of nonminimal couplings, with the same lZi u . Note that, for adiabatic 
initial conditions, the right hand side of (3.18) automatically vanishes and that the matter 
perturbations are in effect adiabatic in all frames. The nonadiabatic pressure term due to 
the change of frame manifests itself only for nonadiabatic initial conditions. 

Let us point out that the equations written in this section include as a special case 
(corresponding to otyi = ot = «b = 0 and ckd,/ = 0) the equations of motion for linear 
perturbations derived in standard models of dark energy (fc-essence) conformally coupled 
with matter (see e.g. [51]). Our results also include the more recent investigations of dis- 
formal couplings between matter, usually CDM, and some standard dark energy (i.e. with 
c^M = = a B = 0) [41—43,47,64]. 

In the general case, eqs. (3.16)-(3.17) can be directly applied to the usual matter 
species, i.e. CDM, baryons, photons and neutrinos and implemented in a numerical code. 
If one wants to study the CMB fluctuations, the fluid approximation is not sufficient for 
photons and neutrinos and must be replaced by a Boltzmann description. Whereas a 
nonminimal coupling of photons is constrained to remain tiny [67], one could envisage a 
nonminimal coupling of neutrinos (see e.g. [68-70]). To deal with this modification, the 
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simplest method would consist in writing the Boltzmann equation in the Jordan frame of 
the neutrinos, where it keeps its usual form. The neutrino-frame gravitational potentials 
appearing in this equation could then be expressed in terms of the gravitational poten¬ 
tials <h and T associated with the baryon-photon frame, by using explicitly the disformal 
transformation between the two frames, as given in Appendix A. 


4 Baryons and coupled CDM 


In this section, we apply the general formalism developed in the previous sections to the 
cosmological era where the dominant matter species are baryons (denoted by the subscript 
b ) and CDM (subscript c). Whereas there exist very stringent constraints on EP violation 
for baryons [26,45], the dark matter sector is much less constrained [19]. For this reason, 
we now assume that the baryons are minimally coupled, i.e. 


ac,fc = 0, «d,() = 0 =+ 76 = 0, (4.1) 

while dark matter is coupled to dark energy via a general metric of the form (2.33). 

For both baryons and CDM, one neglects the pressure and anisotropic stress, so that 
Wb = w c = 0 for the background and 5pb = Sp c = Ob = cr c = 0 for the perturbations 10 . 
The background equations (3.7) and (3.9) take the form 


Pb + 3 Hpb = 0 , (4.2) 

p c + 377(1 - 7 C ) p c = 0 , (4.3) 

PDE + [3(1 + U>de) - «m] -ffpDE = -3i77cPc + HotyiPra • (4.4) 


According to (3.8), the coupling parameter 7 C is related to the CDM conformal and dis¬ 
formal parameters via 


1 

7c = + 


dp, c 

677(1 + «d, c ) ’ 


(4.5) 


4.1 Linear perturbations 

Let us now consider the linear perturbations. The continuity and Euler equations, (3.16) 
and (3.17), reduce to 

k 2 

$b -o Vb = > (4.6) 

cr 

v b = , (4.7) 

k 2 

5 C - kV c = 3(4' + 7 c T77t)' + 2(1 + a D ,c)(ac,c - 37 c )77(4> - n) - a D , c (4 > - vi), (4.8) 

a z 

v c + 3H^f c v c = -4> - 3T77 c 7t . (4.9) 

These equations must be supplemented by the Einstein equations, eqs. (C.2)-(C.5) and 
by the scalar fluctuation equation (C.7). 

It is possible to use a combination of the Einstein equations and of (C.7) to eliminate 
the dependence on n and tt in the above equations in favour of the gravitational potentials. 

10 As shown in Appendix A. 2 , this statement holds in any frame. 


16 




The same procedure has been used in the case of minimally coupled matter in Refs. [17] 
and [14]. In our baryon and coupled CDM system we find a dynamical equation for T of 
the form: 


T + 


Pi fa + faa\ k H tt ^ , PiP* + P 1 P 5 k 2 H + c 2 s a 2 B k A H 2 
Pi + Oi^kjj Pi + ol\ k 2 H 



PiPe,i + Plju\ k 2 H PiP 8 ,i + fa,ia& k H h 
P l + a 2 B k 2 H 01 + Pi + a 2 B k 2 H nVl 


(4.10) 


where ku = k/(aH) and the time-dependent coefficients p a are explicitly given in Ap¬ 
pendix E. They generally differ from those given in Refs. [17] and [14] because the dis- 
formal coupling to dark matter modifies the evolution equation for ir, see eq. (C.7). The 
relation between 4> and T is given by 


a B k z H — T 

«T — «M 


a B 

T 


a 


+ Pi 


4> - T(1 + a T ) ( 1 + a 


. OL T — CKM 


2/3i 


«K — 6as „ 
otBOi H--- Hvi 


(4.11) 


where a = «k + 6ag and we have introduced the time-dependent combination 


£ = o:b( 1 + ay) + ay — olm ■ 


(4.12) 


For ar = a m = 0, this reduces to the familiar relation = T. 

We can also eliminate the dependence on -k from the continuity and Euler equations 
for CDM, eqs. (4.8) and (4.9). For simplicity, we give here the explicit expressions only 
in the case ckm = ckt = 0, for which $ = dt, the generalization being straightforward. In 
this case, the continuity and Euler equations (4.8) and (4.9) become 


S r — 


k 2 P 1 P 2 + £3 k% ■ PiU +Pi^k 2 H + (? s a B ^^-k A H 


a 2 c Pi + a B kjj 


-T + 


Pi + «b k H 


■HV 


3 r> TT \ Pi fa,I + £ 7 ,/ k 2 Plfa,I T £ 9 ,/ k'jj 

V 2 l Pi+a 2 B k 2 H 0l+ Pi + a 2 B k 2 H 


H VI 


v c + 3H-/ c v c = - T - 


37c 

Pi + o^k 2 H 


.f + M , . 2 T 

a —— 1- k H a b'I' 


+ y,¥’ h 


«k - 6«b7/ tt 
otbOi H--- Hvj 


(4.13) 


(4.14) 


where the time-dependent coefficients £ a are given in Appendix E. In this case, where 
= T, eqs. (4.6), (4.7), (4.10), (4.13) and (4.14) form a closed system of equations. 


4.2 Quasi-static approximation 

To investigate late-time cosmology, it is convenient to resort, on sufficiently short scales, 
to the quasi-static limit. This is justified as long as we remain on scales smaller than the 
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sound horizon of dark energy, i.e. k aH/c s [71]. In this limit, the conservation and 
Euler equations for baryons and CDM (eqs. (4.6)-(4.9)) simplify to 


k 2 

- 9^6 = 0 , ( 4 . 15 ) 

a 1 

u 6 = - 4 >, ( 4 . 16 ) 

S c ~~ 2^ = 0, (4.17) 

a z 

v c + 3Hj c v c = - 4 * - 3Hj c tt . ( 4 . 18 ) 


In these equations, all the modifications are encoded in the single parameter j c . Therefore, 
it is not possible to disentangle the conformal and disformal effects. Note that this is due 
to the fact that the nonminimally coupled species is pressureless and that we restrict to 
the quasi-static regime. 

We can then use the generalized Einstein equations to derive the Poisson equation for 
4>. Combining eqs. (C.2) and (C.4) one finds 


-2 4 * — —H 2 Q m { (l + «T + Af ) + [1 + «T + /%(/% + A7)] u c S c } , uii — —— , 

CL Zi 0 Ly\ 


where we have defined the dimensionless parameter 

3\/27 c 


A 7 = 


C S (X 


1/2 


(4.19) 

(4.20) 


which characterizes the strength of the nonminimal coupling of CDM, as well as the 
analogous parameter 




c s ot 1 / 2 


(4.21) 


associated with the modified gravity coefficient £ defined in (4.12). Note that the denom¬ 
inator in the definitions of /3 7 and /3^ is real, since stability requires that c 2 a > 0. 

Moreover, in the quasi-static limit the evolution equation (C.7) for n reduces to a 
constraint equation, which reads 


k 2 

-qTt = 3HQ m 

a z 


P^ujbSb + (/% + (3~()uj c 5 c 

s/^CsOi 1 / 2 


(4.22) 


Substituting (4.19) and (4.22) into the matter equations (4.15)-(4.18), we obtain two 
coupled second-order differential equations for the two density contrasts: 


5b + 2 H5b — —H 2 Q, m {(1 + otT + Af) w f ,5b + [1 + cct + PgiPz + A7)] w c 5 c } , (4.23) 
<5 C + (2 — 3q C )H5 C = —H 2 Q m {[1 + ctT + A^(A^ + A7)] + [l + «t + (A$ + A7) 2 ] 

(4.24) 


In the absence of nonminimal coupling of CDM (Ay = 0), the gravitational coupling of 
both species is modified by the same factor 1 + ckt + Af • I n the absence of modified gravity 


18 





(/?£ = 0 and ay = 0), one finds that the nonminimal coupling of CDM (/3 7 A 0) modifies 
the friction term for 5 C , as well as increases the coefficient in front of 8 C in the second 
equation, whereas all other three coefficients on the right hand sides are unchanged. This 
is the result obtained in the context of coupled dark energy (see e.g. [51]). By contrast, 
if one combines modified gravity (/ 3^ 7 ^ 0) with a nonminimal coupling of CDM, all four 
coefficients on the right hand sides are modified. We leave for the future the detailed study 
of how these new coefficients parametrize the influence of modified gravity on structure 
formation. 

Let us now turn to the two gravitational potentials <f> and T. When considering the 
impact of dark energy on observations, it is often convenient to express the new relations 
between the two potentials 'F and <h and the total matter density fluctuations in terms of 
modifications of the Newton constant. We thus introduce the parameters 

2 M 2 k 2 $ 2 M 2 k 2 ^ 

P < t > =-2- ^ = -2-7— 1 ( 4 - 25 ) 

O PmOm O PmOm 

which are equal to one in the standard case. From eq. (4.19) and an analogous Poisson-like 
equation for T, obtained by combining the Einstein equations, one finds that the above 
parameters are given by 


/i<j) — 1 + ax + /3£(/% + /3 7 cu c 6 c ) , (4.26) 

h'f = 1 + /3b(/% + /3 7 u; c & c ) , (4.27) 


where we have defined 


/3b — 


V2a B 
c s a 1 / 2 


(4.28) 


We have also introduced a time-dependent bias parameter, b c = 5 c /5 m . n 

As the gravitational lensing effect depends on the sum of the two potentials, the rele¬ 
vant quantity parametrizing deviations in weak lensing observables (and equal to two in 
the standard case) is 


gWL — h't + p$ — 2 + ax + (/3 b + /%) (/3^ + /3 7 w c 6 c ). (4.29) 

Thus, the impact of modifications of gravity due to non-vanishing as, «m and ot affects 
observable quantities in the perturbations through ay and the combinations /3 b, /3j. Anal¬ 
ogously, the effect of nonminimal couplings on observations is parameterized by /3 7 only 
(see the next section for the case of a coupled relativistc fluid, where another quantity is 
needed to parameterize the nonminimal coupling). Note that as a consequence of dropping 
time derivatives in the fluctuations of it, the parameter a always appears multiplied by 
(? s . From the definition of the sound speed, eq. (2.39), (? s a is independent of «k, so that 
the latter cannot be constrained by observations in the quasi-static limit [16]. 

When /3 7 = 0, i.e. if CDM is minimally coupled and there are no EP violations, the 
last term inside the parenthesis of eqs. (4.26) and (4.27) drops and these relations simplify 
to 

//vj, = 1 + /3^/3b , /i# = 1 + qt T 0^ ■ (4.30) 

11 In the quasi-static limit, the evolution equations (4.23) and (4.24) are scale independent so that the 
ratios Sb/Sm and 5 c /5 m do not depend on scales. 
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(4.31) 


In this case, the so-called slip parameter becomes (see for instance [14]) 

£ = i + PzPb 

<b 1 + a T + /3| ’ 

By contrast, if there is a non trivial coupling of CDM but gravity itself is not modified, in 
which case we have «b = 0, olm = 0 and ccp = 0 (thus £ = 0), we recover that the Newton 
constant is not modified, /i$ = fly = 1, and that 4> and 4/ are the same as in GR, even if 
CDM is nonminimally coupled, as is the case in usual scenarios of coupled dark energy. In 
general, we find that the situation is much richer when both gravity and matter couplings 
are modified. 


5 Matter and coupled relativistic fluid 


In this section, as another example we consider nonminimally coupled relativistic particles, 
in the fluid approximation. They could represent neutrinos, radiation or warm dark matter 
in the relativistic regime. Baryons and CDM are taken to be minimally coupled, ao,m = 
a c,m = 0. In the Jordan frame of the relativistic fluid, its equation of state parameter 
is given by w T = 1/3 (and c 2 r = 1/3). Thus, in the frame where baryons and CDM are 
minimally coupled, the background and perturbed equations of state are 


w r = 


3(1 + CtD,r) 


5p r = 


Pr 


3(1 + «D,r) 


<5 r + 2a D ,r( < h - 7r) - .. Qp ' r 7T 
1 + «D,r 


(5.1) 


The second relation has been obtained from eq. (3.18), using eq. (2.48) for the sound 
speed. To simplify the treatment, we ignore the anisotropic stress, i.e. er r = 0. 

We are now going to assume that baryons and CDM dominate the gravitational per¬ 
turbations, thus neglecting the backreaction of the relativistic fluid. On small scales, we 
can then resort to the quasi-static approximation. Under these conditions, the evolution 
equations for matter are 

k 2 

K -n«m = 0, h m = -4>. (5.2) 

a z 

For the relativistic fluid, we use eqs. (3.16) and (3.17) and replace 7 r with the expression 
7r = (aD,r — 2i7ao,r)/ [677(1 + «D,r)]- The evolution equations then read 


v T — H 1 


Sr~ 


4 + 3ao, r k 2 
3(1 + «D,r) a- 


;Vr = - 


a D,r 


k 2 


r7T 


(5.3) 


2 + 3«D,r , J. A 

3 -r ~-5 - (5D,r - jD,r) 

1 + 6a D r 


V T + 


4 + 3ctD,r 


3(1 + «D,r) a 2 
- (1 + 2/D,r) < h + 2f D} rTT - 3Hg Dtr 7T . 

(5.4) 


where we have defined 

f = Q D,r 

/D ’ r _ 4 + 3a D , r 


9D,r = /d, r 1 


(1 + 3aD,r)dD,r 
6(1 + aD,r)77aD,r 


(5.5) 


As expected, for ao,r = 0 the effects due to the nonminimal coupling vanish. Note that, 
even if we are in the quasi-static limit, the term in 7 r should be kept in the Euler equation, 
as it is expected to be of the order of Hit and therefore comparable to the other terms. 
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Our assumption that the relativistic fluid does not contribute to the gravitational 
perturbations means that <h and tt are only sourced by CDM and baryons, i.e. 

-^ 2 $ = 7 ;H 2 n m (l + a T + 0i) <5 m , -^ir = 3Htt m ^ 6 m , (5.6) 

a z 2 a z \f2c s a 1 / 2 

which correspond to eqs. (4.19) and (4.22) specialized to the case /3 7 = 0. Therefore, even 
when the extra scalar field is not sourced by matter perturbation (e.g. when j3^ = 0) 
and 7r = 0, the relativistic particles still feel a force F r different from that felt by matter, 
F m , the relative difference being given by (F r — F m )/F m = 2/D, r - This extra force is due 
to the non-adiabatic pressure perturbation 5p nSi d,i in eq. (3.18), induced by the disformal 
coupling out of the Jordan frame of the fluid (see eq. (5.1)). 

To highlight this effect, the Euler and continuity equations can be combined to form 
a second-order differential equation for the density contrast <5 r , sourced by the matter 
perturbations according to eq. (5.6). In the simple case where gravity is not modified, 
i.e. /?£ = clt = 0, we get 


5 r -\- H5 r I 1 + 3aD,r 


! ~ gD,r//D,r 

1 + 3ctD,r 


+ 


k 2 


3a 2 (l + «D,r) 


<5 r = 2 H 2 n r 


1 + CVD,r/4 
1 + CtD,r 


(5.7) 


Unlike in the case of the bayons-CDM fluid, the signature of the disformal coupling here 
is present at the linear level in aD, r , while in eqs. (4.23)-(4.24) it appears at the quadratic 
level through the terms /%/3 7 and f3 2 . 

The main message of this section is that one must define the usual fluid properties 
(such as the equation of state and the speed of sound) in the Jordan frame, where the 
species is minimally coupled to gravity. 


6 Conclusions 

In this work, we have presented an effective description of dark energy and modified grav¬ 
ity, which extends the approach developed in [5] by relaxing the assumption of universal 
coupling of all matter species. Namely, we have allowed each matter species to be asso¬ 
ciated with a specific Jordan frame (or metric), conformally and disformally related to 
the gravitational metric. In this way, we have made connection with a vast sector of the 
literature devoted to the so-called coupled dark energy, with either a conformal coupling 
in most works or a disformal coupling for more recent works. However, in contrast with 
this previous literature, we have considered here a very general description of the gravita¬ 
tional sector, which includes Horndeski’s theories (although not their extensions such as 
G 3 ) instead of general relativity with a quintessence-like scalar field as usually assumed. 

At the level of linear perturbations, the gravitational sector is described by the quadratic 
action given in eq. (2.27), which depends on four time-dependent parameters «k, ub, «m 
and or- As for matter, each species is characterized by two time-dependent parameters, 
ac,i and od j, associated with their conformal and disformal couplings to the gravitational 
metric. This implies that the whole system depends on a total of 2 N$ + 4 time-dependent 
parameters, if Ns species are present. However, there is some arbitrariness in the choice 
of the gravitational metric that is used to define the gravitational and matter sectors. By 
considering a conformal-disformal transformation (2.42) of this metric, the same physical 
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system is characterized by 2N$ + 4 new parameters, which transform according to (2.45) 
and (2.46). Taking into account this “gauge” redundance, which depends on two arbitrary 
parameters, one thus finds that the number of physically relevant parameters is reduced 
to 2 (Ng + 1). 

A very useful result of the present work is the derivation of the linear stability condi¬ 
tions in this very general framework. As the presence of disformal couplings contributes 
to the kinetic energy of the scalar fluctuations, the condition for the absence of ghosts is 
modified. This now requires that a defined in eq. (2.36) is positive. We have checked that 
the stability conditions are invariant under the “gauge transformations” of the parameters 
discussed above. 

We have also written the equations of motion for the linear perturbations and em¬ 
phasized how the usual equations are modified in the presence of modified gravity and 
nonmininal (conformal or disformal) couplings. Special care must be taken when the cho¬ 
sen frame does not coincide with the matter Jordan frame as the relations between matter 
quantities are frame-dependent. For instance, the equation of state parameter, whose 
natural value (e.g. 1/3 for radiation) is defined in the Jordan frame associated with the 
matter species, will be in general different in another frame. 

We have illustrated our formalism by considering two types of scenarios, motivated by 
the already stringent constraints on the nonminimal coupling of ordinary species (baryons 
and photons) to a scalar field. In the first case, we have focused our attention to the 
situation where only CDM is nonminimally coupled to the scalar field. For late cosmology, 
in the quasi-static approximation, we have computed the evolution equations of CDM and 
baryon density contrasts. In the second case, we have assumed that both baryons and 
CDM are minimally coupled but allowing for a relativistic fluid (e.g. neutrinos) with 
nonminimal couplings. These two simple examples illustrate what kind of new effects can 
be produced by the combination of modified gravity and nonminimal couplings. 

It would be interesting to investigate how future observations will be able to constrain 
simultaneously the parameters describing the deviations from GR and those characterizing 
the coupling of matter to this generalized gravitational sector. 
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A Changing frame 

We consider a general disformal transformation of the metric (2.41), which in unitary 
gauge reads 

t dfiu = C{t)g^ v + D(i)<5°<5° , (A.l) 

and study how metric and matter quantities change under this transformation. In terms 
of the two time-dependent parameters C and D, the ADM components of the new metric 
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in unitary gauge are given by 

N 2 = CN 2 - D , N* = N\ h i:j = Chij , (A.2) 

while the intrinsic Ricci scalar and the extrinsic curvature respectively transform as 

N 


R = C~ 1 R , 


K\ = 


N 


k '; + 4 c s ’ 


For the matter sector, the stress-energy tensor in the new frame is 

2 5S t 


rpfAlS _ 

(/) = ’ 


so that 


_ \J ~9 &9a.fi rpOi/3 _ 


N 


-lilts 


{I) V~9 ^9/ms (/) C 5 / 2 VCN 2 - D (/) 


(A.3) 

(A.4) 

(A.5) 


A.l Background 

Let us now set N = 1 and assume a flat FLRW background, ds 2 = — dt 2 + a 2 (t)dx 2 . From 
eq. (A.2), the background metric in the new frame remains flat FLRW, with line element 
ds 2 = — dt 2 + a 2 (t)dx 2 , where we have defined 


t = 


Jfi 


C 


+ OL D 


-dt 


= VCa . 


From this equation, the Hubble rate in the new frame is given by 


~ Ida , 

H = = (1 + ac) 

a dt 


1 + OLY) 


C 


H . 


(A.6) 


(A.7) 


From eq. (A.5), the background energy density and pressure in the two frames are respec¬ 
tively related by 

1 _ \/l + C*D 


PI = 


C 2 VT+^ PI ’ C 2 

where p = —T° 0 and p = T^/3. This implies wj = wj(l + «d)- 

In terms of these quantities the Friedmann equations (2.12) and (2.13) become 

;(Pm + Pde) , (A.9) 


Pi = 


-Pi , 


(A.8) 


H 2 = 


3 M 2 

dH 3 ~ 2 1 

—— -\- —H — --— 

dt 2 2 M 2 


(Pm+PDE) , 


(A.10) 


where M 2 is given by eq. (2.44). Using the expressions above, one can compute the 
relations between the background energy density and pressure of dark energy in the two 
frames. One finds, respectively, 

1 


Pde = 


PDE = 


CVl + 

a/1 + a d 


PDE + 3 M 2 H 2 v/1 + " D 


C 2 


ac(2 + ac) 


ap 

1 + «d 


(A.ll) 


PDE + M 2 ff 2 ' /1 + OD 


(j2 1 (J2 

where as usual a dot denotes a derivative with respect to t. 


( A , a i o H , o"C , «d (1 + ac) 
a c (4 + ac) + 2ac^ + 2 if + H(l + a D ) 


(A.12) 
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A.2 Perturbations 


Let us now study how perturbations transform under disformal transformations. Due to 
the invariance of the gravitational action under disformal transformations, the pertur¬ 
bation equations have the same form in both frames. Thus, we just need to derive the 
relation between perturbation quantities in different frames. 

Introducing n in eq. (A.l) via the time reparametrization (3.1), one finds, up to linear 
order in 7r, 


5oo — C 
90i = C 


900 + T -^-(l + 2ir) 

1 + «D 

901 + -r-- OiTT , 

1 + OD 


9ij = C{1 + 2Ha C 'n-)gij ■ 


2 ac 
1 + o;d 


Hit + 


OiD 

(1 + an) 2 


(A.13) 

(A. 14) 
(A-15) 


Thus, if we start from a perturbed FLRW metric in Newtonian gauge with goi = 0 we end 
up with goi ^ 0 after this transformation. To maintain the Newtonian gauge condition 
goi = 0, we need to supplement the time redefinition (A.6) with a space-dependent shift 
(see Appendix C of [24]), i.e. 

t= [ \ L ^ dt-anJ f 7r(f,x). (A.16) 

J V 1 + «D V 1 + «D 

Then, the new perturbations result from combination of the field redefinition (A.l) and 
this change of coordinates. For the metric in Newtonian gauge, this yields 


= (1 + «d)^ + [ac(l + ocd)H + do] vr , 
T = T - [ac(l + «d) + «d \H it , 

7T = sJC( 1 + ckd) vr • 


For the matter quantities, using eq. (A.5), one finds 
8j = dj + an (4> — -fr) — 

Spi/pi = (1 + an)\spi/pi - wiajj ($ - vr) 

— Wl 


3iLao(l + wi — 7/) + 4Lfac(l + cud) + 


7T . 


3iLao(l T wi 7 /) + 4iLac(l + «d) — -bo - -«d 

2 wj 


TT 


\/C( 1 + «d) r/1 , ^ ! 

VI = 7 —-77—-7 (1 + wi)vi - w/a D vr , 

1 -I- wj{ 1 + an) 

y/1 + Ot£) 

ai= —^ ai - 


(A.17) 

(A.18) 

(A.19) 


(A.20) 


(A.21) 
(A.22) 
(A.23) 


One can relate the pressure and density perturbations via the speed of sound, which is 
defined as the ratio between these two quantities in a coordinate system where the fluid 
is at rest. In the Jordan frame of the fluid, this gives [72] 


5pi 


clidpj - (1 + u)/)/5/ 


o rV/~2 - \ , dwj/dt 

3H(C,J - w,) + T -^ 7 


vr 


(A.24) 
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One can then use eqs. (A.20)-(A.22) to rewrite this equation in a generic frame. This 
yields 


Spi = c 2 sI 5pi - p! 


3H(c s i — wi)( 1 + ol c) + — 


U!j 


+ wi 


(1 + wi)vi + 


«D 


1 + OLY) 


-7T 


- 4Hp I (c 2 sI - wi)a C Ti + Pi(c 2 r + w /) 


2 — 7r) — 


«D 


1 + «D 


-7T 


(A.25) 


where we recall that, from eqs. (A. 8 ) and (2.48), the equation of state parameters and 
sound speeds defined in the two frames are respectively related by wi = (1 + otd)wi and 

cli = (f + anK/. 


B Explicit quadratic action 

B.l Matter action 


For simplicity, we assume that each matter species can be described by a perfect fluid with 
vanishing vorticity (this restriction does not affect the analysis of scalar linear modes). 
It is then easy to write an action in terms of a derivatively coupled scalar field with 
Lagrangian 12 

Sn = 5>, S r = / d 4 x^PPr(Y r ) , Y r = g'^ard^j . (B.l) 

i 


The second-order expansion of the action Si reads 
:( 2 ) 


sp = r ^ dtN 4- ( 


pi 


Sj 


5N 

T 


+ 


1 + (1 + «p,j)wj p! 

(1 + CkD,/) 2 CTI 
1 + (1 + anj)wi pi 
(1 + «d j ) 2 2ct 2 


SN 


S&i f - c 2 jdVh ) + c 2 s I N l di5ai 


x a-2 a“t2„2 

o&j N c s j 9 


(B.2) 


where we have split the scalar field 07 into a background value and its perturbations, 
a 1 = It 1 (t) + 5ai{t,x). The fluid quantities are related to the function T/ (Yj) through 


Pi 




C 2 


\/l + «D ,1 

P'l 

P'j + 2Y t P'/ 


Pi 1 pi — C 2 y 7 ! + an,! ( 2YiPj — Pi) , 


(1 + «D ,l) 


-1 


(B.3) 


where a prime denotes a derivative with respect to the variable Yj. We have omitted in the 
action irrelevant terms that vanish when imposing the background equations of motion. 
For Ci = 1 and «d,/ = 0 we recover the usual expressions for a ^-essence fluid [73,74,76]. 

12 The more general fc-essence type Lagrangian of Refs. [73,74] explicitly depends also on the scalar field. 
Since here we are interested only in the derivative terms, we assume for simplicity that Pi depends only 
on Yi and not on oy. This description implies that each of the fluids is also barotropic [75], i.e. that its 
pressure is a function of its energy density, pi = pi{pi)- 
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B.2 Stability and sound speed of dark energy 

In order to investigate linear stability issues, we need to extract the quadratic action for 
the propagating degrees of freedom. We concentrate on scalar modes as the stability 
conditions of tensors are not modified by the nonminimal coupling of matter. To this end, 
we will expand the total action up to quadratic order in linear scalar fluctuations around a 
FLRW solution and solve the constraints, generalizing the procedure of Refs. [77] and [24], 
The second-order action 

S (2) = g(2) + g(2) ^ ( B . 4 ) 

( 2 ) 

where the gravitational part Sg ; is given in eq. (2.27), governs the dynamics of linear 
scalar fluctuations. Assuming N = 1 without loss of generality, the scalar modes can be 
described in unitary gauge by the metric perturbations [78] 

N = l + 6N, N i = 8 ij d j if>, h ij = a 2 (t)e 2<: 6 ij . (B.5) 

As a consequence, we get 

dVh = 3a 3 C, SK { j = (C - H6N)8) - 5 ik d k djip , (B.6) 

and 

SiRij = -M 2 C - didjC , 8 2 R = ~ [(SC) 2 - 4CS 2 C] • (B.7) 

(The metric perturbations 8N and C and the scalar fluctuation are related to the metric 
perturbations in Newtonian gauge by 5N = — tt, C = — 'h — Hit and if; = a~ 2 n.) 

Substituting these expressions into (B.4), we obtain the second-order action in terms of 
the three scalar quantities 5N, ^ and (. Variation with respect to -0 yields the momentum 
constraint, whose solution reads 


8N 


1 + 1 + ( 1 + ttp,/)w/ nj jgA ; 

1 + «b yH 2 1 + a~D,l cr r J 


(B.8) 


with Qj = pj/(3M 2 H 2 ). 

We do not need the solution of the Hamiltonian constraint, as the longitudinal part 
of the shift t/) only contributes to a boundary term in the action. Replacing the above 
solution into the second-order action and re-expressing the scalar fields perturbations 8aj 
in terms of the gauge invariant variables 


Qi = 6 a - gC , 


(B.9) 


the total second-order action reads, focusing only on the kinetic and spatial gradient parts, 


£( 2 ) = J d 3 xdta 3 ^- 

(jmtjH 


2 r yi | (<9iC) 2 , * 

+gdcd(—— + }^a 


kiH 2 2 2 (0Q/) 2 


b 2 c 2 

l a I c s,I 


Q 2 i-clr 




Qrt-^diQrfiC 


ai 


(B.10) 


26 










with 


5 CC = 

Strit,/ = 

9dCd( 


I 


(1 + «b) 2 

1 


a 


+ ^ -y-(aD,/ - «b) 


I c s,l 


2~(^D,/ — «b) , 


1 + «B "7" C s J 


1 + CtB 


Ha b . . 

777 + 7—;- 1 “ «b (1 + «t) + «T 

H z 1 + CUB 


(B.ll) 

(B.12) 

E /V/ . V 

— (1 + 2«d,7 - «b) , 

(B.13) 


where we have defined the dimensionless coefficients 


a — «k + 6«b “I" 3 y ^ ctD,/ tl/ j kj — 3 


! + (! + aD,/)^/ 


(1 + od,/) 2 


n i 


(B.14) 


Absence of ghosts is ensured by requiring that the matrix of the kinetic coefficients is 
positive definite, which yields the conditions a > 0 and ki > 0. The second condition 
reads pi + (l+a-£>.i)pi > 0, which is the usual Null Energy Condition written in a disformed 
frame. 

Diagonalization of the kinetic-spatial gradient matrix yields the following speed of 
propagation for dark energy, 


a 


c s ~ -i (1 + «b) 


H 

W 


— Qm + Oi T + G!b(1 + ®t) 


cm 3 
+ li + 2 


E[l + (1 + “D ,l)wi 


Absence of gradient instabilities requires c 2 > 0 and c 2 1 > 0. 


(B.15) 


C Perturbation equations 

Here we provide the generalized Einstein equations in the presence of dark energy and 
modifications of gravity. These have been first given in Ref. [5] in terms of the parameters 
of the Effective Field Theory of dark energy [4] and in Refs. [17] (see also [14]) in terms 
of the parameters a a . 


C.l Einstein equations 

Let us defined 

w m = ^2—W! , 7 m = ^—7/, (C.l) 

y' Pm “ Pm 

where 7 / parametrizes the nonminimal coupling of the species /, see definition in eq. (3.8). 
The Hamiltonian constraint ((00) component of the Einstein equation) is 


k 2 


6(1 + aB)L7'I' + (6 — ok + 12ob)H 2 ^ + 2—+ (ok — 6 ob) H 2 tt 


+ 6 


(1 + «b )H + —77 2 O m ( 1 + w m — 7 m ) — 


1 k 2 
3 a 


2 «B 


(C.2) 


Hit = -3n m H 2 5, 


rn 11 5 
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while the momentum constraint ((Oi) components of the Einstein equation) reads 


2'i' + 2(1 + — 2Ha-Q-k + [2 H + 3I7 2 fl m (l + tc m )] n = — 3H 2 Q m (l + tc m )u m . (C.3) 

The traceless part of the ij components of the Einstein equation gives 

$ - (1 + q-t)^ + (cum — (xt)Htt = ~~j^2 ’ (C.4) 

while the trace of the same components gives, using the equation above, 


2T + 2(3 + + 2(1 + ob)H^i 


H — —H O m (l + w m ) + («b H) + (3 + Q!m)(1 + o-b)H~ 


+ 2 

— 2Ho.b 7T + 2 


4* 




H + -H z Q m ( 1 + w m ) - (ag H)' — (3 + ayijotBH 2 


TT 


+ 2 < (3 + cxm)HH + — H 2 Q m [w m — 377(1 + w m — 7 m )] + H > it — 


M 2 


Spn 


2 k 2 \ 

3 a 2 am J ' 

(C.5) 


C.2 Scalar field equation 

The charge Qi is dehned in eq. (3.4). Its perturbation reads 

'dpi 


SQj = 3 H 


7 idi ~ ac,i 


Pi 


— wi6i 


1 + wi k 2 

Pi ~ at d,j 1 , — ~?iP lVl 


1 + «D,/ a 2 


+ 2 H [ac,/(l - 3 wi) — 37 1 ] p/4> - - 
H 


od ,1 


+ «D,/ 


k 2 


$ + 3T — 51 — tt + vji—^tt pi 


+ 


1 + «D,/ 

3 


[ - 2ac,/(l — 3w/)(l + «d,/) + 3w/ao,/ + 37/(2 + an,/)] Pfk 


+ + a -{(“/“D,/ + 7/)^ + [(«D,J + ac,/(l + C£D,l))wi + 7/] H |p/7T . 


(C.6) 


We have checked that this expression agrees with those in the literature (see e.g. [41-43, 
47,64]) in the relevant limits. 13 

The evolution equation for n in the absence of EP violations is given in [5] and can be 
found in [14] in terms of the parameters used in this article. Including the contribution 
of J2i 6 Qi using the above equation, and using the continuity equation, eq. (3.16), this 

13 The expressions for the charge Qi given in eqs. (3.6) and (C.6) are in unitary gauge. To compare to 
those in the literature, one must rescale by a factor tfi, i.e. Qi i-> Qi(f> and 5Qi 5Qi(j) + Qi(4>/4>) it. 
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becomes 


«K 


+ 3 ^2 ckd,/12/ J H~ir + < H 2 (3 + C(m) + H ax + (-ffax) 


3 H 2 ^2 12/ 2ac,/(l - 3«;/)(1 + ckd,/) — 3k;/ ao,/ - 67/(1 + ckd,/) 
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u'2 


Hn 


+ 3< 2H + 3 HH Z 


+ 2HaB 


H ~{3 + ckm) “I - H 


12 m (l + W m ) + ^ w I a D,I^I 
I 

+ 2 H { Ha - o ) +3 H 3 ’^2 [W («d,; + ckc,/(1 + ckd,/)) + 3-07/(1 + k;/ — 7/) 
/ 

k 2 


12/ >7T 


cr 


2# + 3i7"12 m (l + K; m ) + 2i7“ ckb(1 + om) + ckt — ckm + 2 (i7a;j j) 


+ 3i7 2 ^ ' w;/ckd,/12/ >7 r + 6i?aB^ + 77 2 (6ckb — ckk) — 3 H 2 ^ ' ckd,/12/ 

/ / 

+ 3 |^2i7 + 3i7 2 12 m (l + w m ) + 277 2 ckb(3 + ckm) T 2(ckb-H) + 3 H~ ^ ^ wjqidjQj 




+ 16 H + 9i7 2 12 m (l + w m ) + i7 2 (6as — cur) (3 + qm) + 2(9ckb — cxk)H + H(6ctB — ok) 

r 1 'i u2 

— 6 H 2 ^2 [37/ — ac,/(l — 3 wj) (1 + ckd,/)12/ jH<& + 2{ [H (qm — a/)] T — olbH<&} 

+ 9 H 3 "22 | 7/Ar - [«c,/(l + ckd,/) + ckd,/] ~ w i5j\ 112/ = 0 . 

(C.7) 


D Synchronous gauge 

Here we provide the perturbation equations in synchronous gauge, often employed in 
numerical codes, where the perturbed FLRW metric has the form 


ds 2 


—dt 2 + 


l + j/. 


5ij + 


/ kikj 



(h + 67 ) 


dx l dxi . 


(D.l) 


Defining e = a 2 (/i + 6 rj)/k 2 , one can write Newtonian gauge quantities in terms of syn¬ 
chronous gauge ones using the following relations (see for instance [79]), 


$ = e, 

X W r (S 
5p\ = dp) 


T = 7 — He , 
I-Are, 


7r (N) =n (S) + ^ 


jr ( N ) X ( S ) , ■ 

op) =dpj+pie, 


XN) 


= - d f2k 2 -e, 


(D.2) 


where we have introduced the divergence of the velocity, 6j = —k 2 vj/a. (The anisotropic 
stress is gauge invariant.) We can then use the above relations to rewrite eqs. (C.2)-(C.5) 
in synchronous gauge. To do this, we use conformal time, r = J dt/a, and denote by a 
prime the derivative with respect to it. Rescaling the scalar fluctuation n by the conformal 
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factor, 7 r —>• 7r/a, and defining the conformal Hubble rate as H. = a'/ 0 ; one obtains 
((00) component) 

2/c 2 ry - H{ 1 + ae)^ - % 2 (6aB - ok)^ + 9% 2 fi m (l + u; m - 7 m ) 

1 a 2 (D ' 3) 

— 'H 2 (6 — «k + 12«b) + 6ft'(l + ob) — 2A 2 as Hn = — > 


((Of) component) 


2r/ - 2^a B vr' + ^ - 2(1 + a B )H 2 + 3ft 2 f2 m (l + u; m )J tt = w (p m + p m )-^ , (D.4) 


((ij)-traceless) 


2k 2 


h" + 6??" + H{ 2 + otyfjih' + 6r/) — 2 Ai 2 (1 + a T )?? - 2A’ 2 "H (cut - «m)7t = ^p a ^ > ( D - 5 ) 
and ((ij)-trace) 

h" + H(2 + aM)h! — 2k 2 (l + «t)?7 + GaB'Hn" 
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6% 2 a B (3 + a M ) + 6(a B ^)' - 9H 2 H m (l + w m ) - 6(ft' - TA 2 ) 


7T 


+ |6"H“[2 + aM + «b(2 + qm)] + 6(as — + 6(aB - H) / — 2 A' 2 (qt — qm) 

l 1 a 2 

(1 - 3w m )(l + w m ) + 3w m j m + w' m /U - = -3jj7jSp m . 


(D.6) 


- 9n 2 n r 


In synchronous gauge, the evolution equation for the scalar fluctuation, eq. (C.7), reads 
T~L 2 ^ax + 3 a D,/^/^ 7r,/ + |% 2 aK(2 + qm) + H'chk + (olyJ-L)'+ 


3 H 2 ^ f1/ 2ac,/(l — 3«;/)(1 + ao,/) — ao,/(l + 3 wi) — 67/(1 + ao,/) 


Un 


— 2k 2 < "H 2 [aBdM + «t — «m — l] + (aB - H) / + %' + — 1 + rc/(l + ao,/) 


H/ >7T 


+ — 6aB«M + «k(1 + ckm)] + 377 2 'H / [—4 + ax — 2aB(3 — qm)] 

+ 6(1 + a B )U' 2 - 7A 3 (6a' B - a' K ) + 6 U{a B U')' 

+ 3"H 2 fl/ — 27A 2 ac,/(l — 3 k;/)( 1 + a^,/) + 'H , ^3 + an,/ + 3 k;/( 1 + ao,/)j 
1 L 

+ 3 ? 7 2 ^ 7 /(l + 2 ao,/) — (1 — 37 /)(l + k;/ — 7 /)^ + 37Aw/^aD,/ + ac,/(l + od,/)) 

— HaBh" — |'H / + "H 2 ob(1 + qm) — 1 + {ocbT-L)' + 7 ^ 7 A 2 1 + wj( 1 + ao,/) Hr |h/ 
+ 2A’ 2 "H(aM — ot)7 = 0. 

(D.7) 
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The continuity and Euler equations for matter become, respectively, 

d'j + 3%(1 + ac,/)(l + cud,/) ~ + (1 + wj)9i — anjir" 

-H[ — 2a c ,/(l — 3u;/)(l + «d,/) + «d,/(1 + 3 wj + 67 /) + 97 /] 7 r' 
+ | — 2"H 2 (1 + Q!d,/)( 37 / — ac,/(l — 3 wi)) + wia^jk 2 
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1 r 
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+ - 1 + wi(l + a D j) h' = 0, 
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e, - , + 
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3^7/ 2 

-K 7T . 


%(l + iuj)J p/(l + w;/) 3a 2 ///(l + wi) (1 + wi) 
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E Definitions of the parameters 

The coefficients f3 a appearing in eqs. (4.10) and (4.11) are defined as 

. _ 3 0 1 JH \ 9 
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(E.5) 



(E.6) 



(E.7) 



(E.8) 



(E.9) 


@9,1 — — (4 + 3 c 2 + «m + Q!t) + ^3 j 

where we remind that £ = q:b(1 + cut) + — ckm, a = ctK + 6a| + 3aD iC n c and we have 

defined a = «k + 6 a| = a — 3aD, c ^c- Setting ao, c = otc, c = 7 c = 0 in these equations, 
one recovers the expressions of [17] and [14]. 
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For ctm = ckt = 0, the coefficients 4a appearing in eqs. (4.13) and (4.14) are defined as 
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